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We consider the Peierls instability due to the interaction 
of electrons with both acoustic and optical phonons. We sug- 
gest that such a transition takes place in carbon nanotubes 
with small radius. The topological excitations and the tem- 
perature dependence of the conductivity resulting from the 
electron-lattice interactions are considered. 

PACS numbers: 63.20.Kr, 63.70.+h, 05.45.Yv, 72.80.Rj 

It is well-known that half-filled conducting chains are 
unstable against a doubling of the unit cell, where the 
bond length between neighboring lattice sites alternates 
along the chain (dimerization) . The resulting alternation 
of the electron hopping amplitudes opens a gap in the 
electron spectrum, turning the system into a semiconduc- 
tor (Peierls instability) . The phonons relevant for this 
instability have wave vector q ^ T:/a (ais the distance be- 
tween neigboring lattice sites). These phonons backscat- 
ter electrons from the left part of the one-dimensional 
Fermi surface (wave vector k ~ ~5^) to the right one 
(fc ~ ^) and vice versa. It is usually assumed that above 
the phase transition temperature Tc these phonons have 
a finite frequency and we will, therefore, refer to them as 
optical phonons. Due to the mixing with the low-energy 
electron-hole excitations, these phonons soften and at Tc 
their frequency vanishes (giant Kohn anomaly) |Q. 

In this Letter we show that acoustic phonons (i.e. 
phonons of small wave vector q and small frequency 
<^a('?) = '^'ol'zl) may lead to a similar instability and that 
in the presence of both optical and acoustic phonons, the 
acoustic ones are the first to soften, independent of the 
coupling strengths involved. The acoustic component im- 
mobilizes the topological excitations (kinks) over the or- 
dered state, while the softening leads to a sharp increase 
of the electrical resistivity just above Tc- 

At first sight, acoustic phonons cannot open an elec- 
tronic gap, as they only have small wave vectors, which 
does not allow for backscattering. This argument no 
longer holds, however, in lattices where the unit cell al- 
ready contains two sites above Tc (where the electron 
hopping amplitude is still uniform along the chain), as 
in that situation the phonon wave vector is only con- 
served up to a multiple of For instance, due to 
its zigzag structure the undistorted frans-polyacetylene 
chain has lattice period 2a (Fig. 1(a)). In Fig. 1(b) it 
is shown how a small-g distortion perpendicular to the 
chain leads to bond-length alternation, which results in 
backscattering of electrons. We note that in the strictly 
one-dimensional Su-Schrieffer-Heeger model ||] of trans- 



poly acetylene, this acoustic backscattering cannot occur. 

Another example is provided by an armchair single- 
wall carbon nanotube (CNT) (Fig. 1(c)). The undis- 
torted armchair CNT has two half-filled one-dimensional 
electronic bands and is metallic [Q. Its unit cell con- 
tains two inequivalent carbon atoms, which gives rise to 
backscattering of electrons off acoustic phonons. The rel- 
evant acoustic phonons are the so-called twistons. It is 
easy to see that a uniform twist of the armchair CNT re- 
sults in an alternation of the carbon-carbon bond lengths 
along the nanontubc (Fig. 1(d)), which opens a gap in 
the spectrum of electron excitations. The scattering of 
electrons on twistons has recently been used to explain 
the linear temperature dependence of the resisitivity of a 
single-wall CNT g. 
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FIG. 1. The zigzag structure of frans-polyacetylene with 
two carbon atoms per unit cell, without lattice distortion (a) 
and with a small-wave-vector (acoustic) deformation perpen- 
dicular to the chain (b). Similarly, the undistorted armchair 
CNT (c) consists of connected zigzag chains in the tube di- 
rection. A twist deformation of the CNT leads to an out of 
phase alternation of the bond lengths in neighbouring chains. 

The Hamiltonian of our model has the form, 



^ dx 



(1) 



where the first term describes the kinetic energy of elec- 
trons close to the Fermi points and the backscattering 
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of the electrons due to the lattice modulation induced 
by both optical and acoustic phonons. Here, (Ti and (T3 
are the Pauli matrices and we set h = I. The spinor 
^ncr = (y!^"") describes electrons moving to the right/left 
with the Fermi velocity vp. We have omitted scattering 
contributions that leave the electron on the same side of 
the Fermi surface {R ^ R and L ^ L), as these processes 
do not lead to an instability. The index n — 1, . . . ,Nh de- 
notes the electron band. The second and the third terms 
in Eq.(^ describe the energy of the optical (i = o) and 
acoustic {i — a) phonons, 



dxAj 
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dx 7r|, 



(2) 



TrXiVp 

where the dimcnsionless coupling constant absorbs a 
factor Nh for the number of electron bands, and pi is the 
mass density. The momentum density tt^ is related to the 
phonon amplitude Ui by: TTi{x,t) — pidui{x,t)/dt. 

In the continuum description of the optical and acous- 
tic phonons coupled to electrons, is related to the 
phonon amplitude Ui by 



Ao oc Uo{x,t) for optical phonons, 
Aq oc -^Ua{x,t) for acoustic phonons. 



(3) 



The difference between the optical and acoustic phonons 
is due to the fact that the "order parameter" A incor- 
porates also the electron-phonon coupling constant, as is 
clear from Eq. (|]). While the couphng ga{q) to the acous- 
tic phonon with wave vector q is proportional to g, the 
coupling to the optical phonon, goi,q), is approximately 
constant. At the same time, Eq. (|^) gives the linear dis- 
persion for acoustic phonons, LUa — vo\q\, and a finite 
frequency ujo for optical phonons. We note that, though 
the coupling to acoustic phonons ga{q) is small for small 
q, the actual strength of the interaction is given by the 
dimcnsionless coupling Ai oc gf /{vpuj'^), which is finite 
both for optical and acoustic phonons. 

For Aa = and Nb = 1 the Hamiltonian Eq. (0) 
coincides with the Hamiltonian of the TLM model g|, 
which is the continuum version of the Su-Schrieffer- 
Heeger model Q of irans-polyacetylene with one electron 
band. On the other hand, for Aq — and Nb = 2, the 
Hamiltonian Eq. (|^) is equivalent to the Hamiltonian de- 
scribing the electron-twiston interactions in the armchair 
CNTs 01. 

Phase transition. — First we consider the temperature 
dependence of the optical and acoustic phonon frequen- 
cies for T > Tc- The bare phonon propagation Dq and the 
the coupling V to electrons is formally written in terms 
of 2 X 2-matrices, 



Dn = 



Dl 
D^^ 



and 



9 a 9a 9a 9o 
9o9*a 9o9*o 



where Do contains the bare acoustic and optical phonon 
propagators D'^ and D^. The coupling V of the phonons 



to electrons has two effects: It mixes the bare optical and 
acoustic phonon modes and renormalizes their frequen- 
cies due to an electron-hole excitation that is described 
by the bare vacuum polarization P = P(q,w,T) ||]. 
Within the random phase approximation, the propaga- 
tion of dressed phonons is described by a matrix D that 
solves the Dyson equation D = Dq + DqVD. The 
renormalized phonon frequencies (Do (9) and uJaiq) are 
found from the poles of det{D), or equivalently by solv- 
ing det{D^^ - V) = 0. Using that a)a(g) < T and 
'^aiq) ^ 'Do(g), we eventually find: 



^liq) ^^l{q) 



1 ~ (Ao + Ag) H-fW/nT) 
1 - Ao ln(7iy/7rT) 



(4) 



for the acoustic phonon and, assuming that 1^0(9) ^ T , 

Col{q)=u,l{q){l~\oH-iW/iTT)) (5) 

for the optical phonon. Here, W is the energy cut-off of 
the order of the electron band width and 7 = 1.781072... 
is Euler's constant. 

We thus see, that the expression for the renormalized 
optical frequency is independent of A^ and is the same as 
in the absence of the coupling to acoustic phonons. On 
the other hand, the renormalized acoustic phonon fre- 
quency depends on the sum of Aq and Aq. As a result, 
the acoustic phonons "soften" first, at the critical tem- 
perature given by 



TT V Aa + Ac 



(6) 



As ojaiq = 0) = at all temperatures, the "softening" 
in this case means vanishing of the acoustic phonon ve- 
locity at T = Tc. Thus, no matter how much the opti- 
cal coupling constant is larger than the acoustic coupling 
constant, it is always the velocity of the acoustic phonon 
that becomes zero at the transition temperature, whereas 
the optical phonon frequency stays finite at T = Tg. This 
is a consequence of the mixing of the optical and acoustic 
phonons due to their interactions with electrons, which 
results in a repulsion between the frequencies of the two 
modes. As a result, the optical and acoustic branches 
can never cross and the singularity at Tc always occurs 
in the lower, i.e. acoustic, branch. A similar effect takes 
place in some ferroelectrics, in which the sound veloc- 
ity vanishes at the transition temperature because of the 
mixing of the soft mode, describing the ferroelectric dis- 
placement of ions, to acoustic phonons We also note 
that, at first sight, Eq. (||) resembles the result for the 
Peierls temperature obtained in Ref. ||l^. It should be 
kept in mind, however, that in that paper the various 
coupling constants A^ correspond to contributions from 
scattering within different electron bands and are not as- 
sociated with the presence of several phonon modes. The 
additive effect of the number of electron bands is implicit 
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in our result through the fact that both Aq and are 
proportional to Ni,. 

Next we discuss the ordered state below Tc. Because 
the lattice distortions, corresponding to the optical and 
acoustic phonons, are coupled due to electron-phonon 
interactions, both order parameters, Aq(x) and Ao(a;), 
should appear below the transition temperature. In the 
mean field treatment of the lattice one has to minimize 
the total free energy of the model Eq. (|l|) with respect to 
the two order parameters. The solutions of the resulting 
self-consistency equations have the following properties: 
(i) The sum of the optical and acoustic order parameters, 
A(a;) = Ao(cc) + Aa{x), satisfies the same Bogoliubov-de 
Gennes equations as the order parameter of the TLM 
model 1^ with a single phonon mode, but with the cou- 
pling constant A = Ao -|- Aa; (ii) The optical (i = o) and 
acoustic (i — a) order parameters are proportional to 
A{x): 

A,{x) = ^A(x). 

Thus, the solutions of the TLM model can be used to 
study the ordered state in our model. In particular, the 
order parameter for the homogeneous solution A{x) = Aq 
at zero temperature has the value Aq = -Tc, with Tc as 
inEq. §). 

The homogeneous optical order parameter, Aq, is 
a "frozen" phonon mode, corresponding to a uniform 
dimerization of the chain. On the other hand, for acous- 
tic phonons, the constant Aa corresponds to the linearly 
growing amplitude of ionic displacements: Ua{x) oc AaX. 
Such a lattice distortion is not a "frozen" phonon mode, 
as it corresponds to large deviations of ions from their 
equilibrium positions in the high-temperature phase. In 
the case of CNTs, Aq — const, corresponds to a uniform 
twist of the CNT (see Fig. |l|). Such a twist was, in fact, 
recently observed using STM g. 

Topological Excitations. — In the SSH model, where 
electrons interact with one optical phonon mode, domain 
walls (kinks) in the order parameter, corresponding to a 
change of sign of the lattice dimerization along the chain, 
constitute an interesting class of excitations In the 
continuum model, the analytical expression for the kink 
is given by |^: 

A(x) = Aotanh-^, (7) 

where fn = ^ is the correlation length. 

The kink Eq. ([7|) also is a consistent solution of the 
Bogoliubov-de Gennes equations for the total order pa- 
rameter in our model with both optical and acoustic 
phonons. For this kink, the optical and acoustic lattice 
distortions take the form: 

Uo{x) = -Uotanh^, 

(8) 

Ua{x) — Ua In cosh ^. 



Near the kink both the lattice dimerzation, described by 
Uo{x), and the derivative of Ua{x) change sign (see Fig. ||). 
In the case of a CNT, the latter corresponds to the change 
of the sign of the twist-angle from — (ua/Co) to -l-(Ma/^o)- 




FIG. 2. Schematic plot of the optical {uo{x)) and acoustic 
(ua{x)) lattice distortions Eq. (^) corresponding to the topo- 
logical excitation Eq. (j^). 

The creation energy of the kink in the model Eq.dl[) 
is the same as for the kink in the TLM model M: 
Ek = Ni,{2Ao/tt). The dynamical properties of the kinks 
in these two models are, however, quite different. In 
the TLM model the kink can propagate with velocity 
V along the chain, without changing its profile. This re- 
sults from the independence of the kink energy of its po- 
sition and the fact that the kinetic energy density of the 

moving kink, Uo{x,t) = Uq tanh decays exponen- 

tially at distances larger than from the kink. Thus, 
the mass of the kink (defined by equating the kinetic en- 
ergy of the moving kink to iAfw^) is finite and for trans- 
polyacetylene it was estimated to be ~ 6me On the 
other hand, the motion of the kink in the model with 
acoustic phonons would result in a constant kinetic en- 
ergy at distances larger than from the kink, as follows 
from the substitution of Ua(x) in Eq.(||) by Ua{x — vt). 
Thus, the mass of such a kink is proportional to the chain 
length L. This relates to the fact that a translation of the 
kink configuration Eq.(||) changes the coordinates of the 
chain ends, so that a shift of the kink induces a motion 
of the entire chain. Thus, in our model, isolated kinks 
cannot propagate along the chain. Kinks can propagate 
without affecting the chain ends only together with an- 
tikinks, the corresponding time-dependent lattice config- 
uration being Unix, t) — Ua{x — vt) — Ua{x + R — vt) [cf. 
Eq. (^)]. The mass of such a pair is proportional to the 
distance R between the kink and antikink. 

Electrical conductivity. — Next we consider the tem- 
perature dependence of the electrical conductivity, a{T), 
above Tc. For T Wq, the contribution of the opti- 
cal phonons to the conductivity can be neglected. For 
Aa <C Ao this inequality may not be fulfilled close to 
Tc, as the renormalizcd optical phonon frequency at Tc, 
ujo — Wo-\/Aa/(Aa + Ao), may then be small. However, as 
we argue below, for CNTs, Aa ~ Ao, in which case there 
is no dramatic softening of the optical phonon and the 
conductivity close to Tc is also dominated by the electron 
backscattcring off acoustic phonons. 
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The electrical conductivity is given by <t{T) — 
{—Ae^vp /-k) J dkTk{dnF{vFk)/dk), where np is the 
Fermi distribution and is the transport lifetime, which 
depends on the electron wave vector k. Using Fermi's 
Golden Rule and accounting for the renormalization 
Eq. fl) of the acoustic phonon frequency, we obtain 
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FIG. 3. The electrical resistivity p — l/cr as a, function 
of temperature. The electron-phonon coupling strengths are 
chosen to be Xa = Xo = 0.05. 

Fig. ^ shows the temperature dependence of the elec- 
trical resistivity, p{T) = l/cr(r), calculated for T > 
using Eq.(||). At temperatures T Tc, the resistivity 
decreases linearly with temperature, p(T) oc XaT, as was 
also found in Refs. djlj]. However, close to Tc, due to the 
vanishing of the acoustic phonon velocity at the critical 
temperature, the resistivity strongly increases up to some 
finite value, p{Tc). This behavior is very similar to the 
one observed for bundles of single- wall CNTs where 
the cross-over between the linear decrease and the sharp 
upturn of the resistivity occurs at T* 10 — 100 K. 

Usually, the Peierls instability is thought to be irrel- 
evant for CNTs, because the electron-phonon coupling 
constant is inversely proportional to the number of zigzag 
chains, 27V, in the armchair CNT, making Tc negligibly 
small for iV > 4 [Q. We note, however, that previously 
only the interaction with optical phonons, e.g., the de- 
formations resulting in a Kekule or quinoid type struc- 
ture of the CNT lattice, was considered When we 
take into account the scattering of electrons on acoustic 
phonons (twistons), we obtain a higher value of the cou- 
pling constant A (= Ao -I- Xa), which gives a much higher 
value of Tc. In particular, using typical parameters for 
graphite ||], we obtain Xa — 0.23/N. Furthermore, if one 
assumes that the electron hopping amplitudes and the 
spring constants only depend on the carbon-carbon bond 
length, one finds Ao ~ Ao. For CNTs with = 4, which 
were recently discussed in the context of possible techno- 
logical applications |l^ , the critical temperature given by 



Eq. (||) is then ~ 10 K, instead of only a few mK for the 
Peierls transition temperature without acoustic phonons. 
We note, however, that the interaction of the CNTs with 
their environement (other CNTs and the substrate) may, 
in principle, suppress the Peierls instability. 

The values of the electron-phonon coupling constants 
for a CNT could, in principle, be obtained from the form 
of the resistivity curve {Xa determines the slope of p{T) 
at high temperatures, whereas Tc depends on Xa -\- Xo). 
However, so far, such measurements have only been per- 
formed on bundles of nanotubes, in which one cannot 
control the radius and chirality of the CNTs. Therefore, 
measuring the resistivity of a single CNT would be very 
desirable. 

Conclusions. — In this Letter we studied the Peierls 
transition resulting from electrons interacting with both 
optical and acoustic phonons. We found that, inde- 
pendent of the electron-phonon coupling constants, the 
acoustic phonon velocity vanishes at the critical temper- 
ature, whereas the optical phonon frequency remains fi- 
nite. This is different from the conventional Peierls sce- 
nario, in which the optical phonon softens at Tc. We sug- 
gested that such a transition may take place in CNTs of 
small radius, leading to a static twist below Tc. The tem- 
perature dependence of the resistivity calculated within 
our model qualitatively agrees with experimental data 
on CNTs. Unlike the topological excitations in trans- 
polyacetylene, the kinks in our model are immobile even 
in isolated chains and can propagate only in pairs. Fi- 
nally we note that, though we used CNTs as a specific 
example, we believe that our results have a more general 
significance. 
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